Collapse arrest and soliton stabilization in nonlocal nonlinear media 
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We investigate the properties of localized waves in systems governed by nonlocal nonlinear Schrodinger 
type equations. We prove rigorously by bounding the Hamiltonian that nonlocality of the nonlinearity 
prevents collapse in, e.g., Bose- Einstein condensates and optical Kerr media in all physical dimensions. 
The nonlocal nonlinear response must be symmetric, but can be of completely arbitrary shape. We use 
variational techniques to find the soliton solutions and illustrate the stabilizing effect of nonlocality. 



PACS numbers: 42.65.Jx, 42.65. Tg, 42.65.Sf 



Collapse is a fundamental physical phenomenon well- 
known in the theory of waves in nonlinear media. It refers 
to the situation when strong contraction or self-focusing 
of a wave leads to a catastrophic increase or blow-up of 
its amplitude after a finite time or propagation distance 
(see for reviews). Wave-collapse has been observed 
in plasma waves pj (the famous Langmuir wave-collapse) , 
electromagnetic waves or laser beams [|| (also called self- 
focusing), Bose-Einstein Condensates (BEC's) or matter 
waves |J , and even capillary-gravity waves on deep water 
0. The effect of collapse appears also in astrophysics, 
where the gravitational attraction plays the same role 
as the self-focusing of electromagnetic waves, tending to 
compress stars of sufficient mass, eventually leading to 
their collapse into a black hole Q . 

Typically the contraction must be able to act freely in 
two or more dimensions to be strong enough to generate 
a collapse. Moreover, the so-called norm, which is the 
power for electromagnetic and plasma waves, the atom 
density for BEC's, and the mass for stars, must be above 
a certain critical value for a collapse to occur. Most com- 
monly the collapse has been discussed in the context of 
the nonlinear Schrodinger (NLS) equation, which is a uni- 
versal model for dispersive (or diffractive) weakly nonlin- 
ear physical systems S. The NLS equation models, e.g., 
all systems mentioned above, in which a wave-collapse 
has been predicted and verified experimentally. 

The collapse singularity is an artifact of the model and 
signals the limit of its validity. Close to the singularity, 
when the amplitude is extremely high and the temporal- 
and spatial scales are extremely short, new physical pro- 
cesses will come into play • A common effect is non- 
linear dissipation, such as two-photon absorption of elec- 
tromagnetic waves and inelastic two- and three-body re- 
combination for matter waves, which efficiently absorbs 
the collapsing part of the wave. Thus collapse acts as 
an effective nonlinear loss mechanism, as is well-known 
in, e.g., Langmuir turbulence Jl(| and BEC's pH|. Ef- 



fects, such as discreteness, saturation of the nonlinearity, 
and nonparaxiality, will also eliminate the possibility of 
a collapse (see Q). In contrast, effects such as weak 
linear loss [pf, t emperature fluctuations |^2| , and spatial 
incoherence [la ], cannot eliminate the collapse. In any 
case the collapse effect represents a strong mechanism 
for energy localization, which is important to study to 
understand the properties of a given physical system. 

The inherent nonlocal character of the nonlinearity has 
attracted considerable interest as a means of eliminating 
collapse and stabilizing multidimensional solitary waves. 
Nonlocality appears naturally in optical systems with a 
thermal [jl4| or diffusive jl5| type of nonlinearity. Nonlo- 
cality is also known to influence thepropagation of elec- 
tromagnetic waves in plasmas [fl"^— [ig|] and plays an im- 
portant role in the theory of BEC's, where it accounts 
for the finite-range many-body interaction IllLptj p2| . 



In this work we consider NLS equations with a gen- 
eral nonlocal form of the nonlinearity. Turitsyn proved 
the absence of collapse for three particular shapes of the 
nonlocal nonlinear response |^] . The analysis of the col- 
lapse conditions for general response functions is difficult 
and has been carried out only numerically However, 
in many systems, such as BEC's, one has no knowledge 
of the particular response function, and thus it is impor- 
tant to maintain its generality in the model. Here we 
prove rigorously that nonlocality eliminates collapse in 
all physical dimensions for arbitrary shapes of the nonlo- 
cal response, as long as it is symmetric. 

We consider the evolution of a wave field u = u{r) = 
u(f,r) described by the general nonlocal NLS equation 

du 

i— + V 2 u-V(r)u + N(I)u = 0, (1) 

OT 

where V=V (r ) is an external (linear) confining potential, 
I=I(r)~I(r,T)—\u\ 2 , t is the evolution coordinate and 
f=(ri, r 2 , rz) spans a D-dimensional "transverse" coordi- 
nate space. The nonlinear term N=N(I) is represented 
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in the general nonlocal form 



N(I) = / R(r' - r)I(f')df", 



(2) 



where we assume the response function R(r) to be real, 
localized, nonsingular, and symmetric. For localized 
waves Eq. ([!]) conserve the power (in optics) or number 
of atoms (for BEC) P and the Hamiltonian H, 



P 



Idr, H 



\Vu\ 2 + VI -^-}dr. (3) 



In optics u is the envelope of the electric field with in- 
tensity / and V represents a guiding structure (waveg- 
uide). Here Eq. (^) represents a general phenomenologi- 
cal model for self-focusing Kerr-like media, in which the 
change in the refractive index induced by an optical beam 
involves a transport process. This includes heat conduc- 
tion in materials with a thermal nonlinearity ]T^ | or dif- 
fusion of molecules or atoms in atomic vapours jl5| . A 
nonlocal response in the form (^J) appears naturally due 
to many body interaction processes in the description of 
BEC's Opj,B2,B4| , if the assumption of a zero-range 



For BEC's with a 
) is the nonlocal Gross- 



interaction-potcntial is relaxed 
negative scattering length Eq. (|l 
Pitaevskii (G-P) equation for the collective wave function 
u (t is time), with / representing the density of atoms 
and V representing the magnetic trap. 

In the limit when the response function is a delta- 
function, R(r)=5(\r\), the nonlinear response is local (see 
Fig. la) and simply given by 



N(I) = I, 



(4) 



as in local optical Kerr media described by the standard 
NLS equation and in BEC's described by the standard G- 
P equation. In this local limit multidimensional optical 
beams with a power higher than a certain critical value 
(P > 11.69) will experience unbounded self- focusing and 
collapse after a finite propagation distance It is 

also well-known that BEC's would collapse when the to- 
tal number of atoms is larger than a critical number p| . 

With increasing width of the response function R(f) 
the wave intensity in the vicinity of the point r also con- 
tributes to the nonlinear response at that point. In case 
of weak nonlocality, when R(r) is much narrower than the 
width of the beam (see Fig. lb), one can expand I(r') 
around r'—f and obtain the simplified model 



N(I) = I + 7 V 2 /, 



1 

7= 2 



2 R(r)df 



(5) 



where the small positive definite parameter 7 measures 
the relative width of the nonlocal response. The diffusion 
type model (|^) of the nonlocal nonlinearity is a model in 
its own right in plasma physics, where 7 can take any 
sign fl^ Ji7| . I t was also applied to BEC's p2fl , nonlin- 
ear optics }25| , and energy transfer in biomolecules [ p6[ . 
In weakly nonlocal media with N(I) = I + 7V 2 / it is 



straightforward to show that collapse cannot occur. This 
was first done for plasmas |l^| , and later for BEC's [p2[ . 
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FIG. 1. Degrees of nonlocality, as given by the relative 
width of the response function R and the intensity profile 
/ in the x-plane. Shown is the local (a), the weakly nonlocal 
(b), the general (c), and the strongly nonlocal (d) response. 

In the limit of a strongly nonlocal response much 
broader than the characteristic width of the wave func- 
tion (see Fig. Id), one can instead expand the response 
function and obtain (to lowest order) 



N(I)^P(R Q + R 2 r 2 ). 



(6) 



where R a = i?(0) and R 2 = ±V 2 R(6). The evolution of 
optical beams in such a strongly nonlocal medium was 
considered in |27j . Since this relation is linear, the highly 
nonlinear effect of collapse cannot occur. 

So in the two extreme limits of a weakly and highly 
nonlocal nonlinear response the collapse is prevented. 
For arbitrary degree of nonlocality it is difficult to prove 
anything rigorously. Just saying that the dynamics is 
described by either the weakly nonlocal model (||) or the 
linear oscillator model @, which both have no collapse, 
is not enough. As it is well-known from studies of gen- 
eral NLS equations, the typical singularity is a so-called 
blow-up featuring the amplitude locally going to infinity 
on a broad background localized structure [Q. Such a 
two-scale field distribution, which was also recently ob- 
served in BEC's ^lj, is clearly described by neither of 
the two simple limiting systems. 

The stabilizing effect of nonlocality of an arbitrary de- 
gree was proven by Turitsyn for three specific examples, 
including Coulomb interaction (R(f)=l/\r\) He 
bounded the Hamiltonian from below for fixed power, 
which proves that a collapse cannot occur and that the 
soliton solutions are stable in a Lyapunov sence. Here 
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we consider the general case of arbitrarily shaped, non- 
singular response functions and prove rigorously that the 
Hamiltonian is bounded from below in all dimensions. 

Introducing the D-dimcnsional Fourier transform (de- 
noted with a tilde) and its inverse 



I(k) 



I{f)e lkr df, 1(f) 



1 



(27T) 



I(k)e- ik - f dk, 
(7) 



it is straightforward to show that for N(I) given by 
Eq. (0) the following relations hold 



I(k)\ = 
NIdf 



I(r)e lkr df 
1 



(27T) 



D 



< j Idf=P, (8) 
R(k)\I(k)\ 2 dk. (9) 



For response functions with a positive definite spectrum, 
R(k ) > 0, the Hamiltonian is thus bounded as follows 



H > 



\Vu\\ 2 2 - l -R P 2 



(10) 



where ||u||2 = / |u| p df and we have used the fact that 
the minimum of V(f) is zero. 

The inequality (|l(]) is the main result of this letter. It 
shows that for all symmetric response functions with a 
positive definite Fourier spectrum and a finite value at 
the center, the Hamiltonian is bounded from below by 
the conserved quantity -t^RqP 2 , or conversely, that the 
gradient norm ||Vu||| is bounded from above by the con- 
served quantity H + ^RoP 2 . This represents a rigorous 
proof that a collapse with the wave- amplitude locally go- 
ing to infinity cannot occur in EEC's, plasma, or optical 
Kerr media with a nonlocal nonlinear response, for any 
physically reasonable response function. 

The stabilising effect of the nonlocality can be further 
illustrated by the properties of the stationary solutions 
of Eqs. (0). As a simple example we consider nonlocal 
optical bulk Kerr media with a Gaussian response 



1 



exp 



(11) 



The ground-state stationary solutions are then radially 
symmetric bell-shaped, nodeless solutions of the form 
u(f,z) — <f>(r) exp(zAz), where the profile 4>(r) is found 
from the Euler-Lagrange equations for the Lagrangian 



L = 



M 2 + |V# 



l 



dr. 



(12) 



To capture the main physical effects we use the approx- 
imate variational technique with a Gaussian trial profile 
4>(r) = aexp[— (r//3) 2 ], in view of the fact that the Gaus- 
sian profile is an exact solution in the strongly nonlocal 
limit with N(I) given by the parabolic potential (^). 



Inserting this ansatz into the Lagrangian (|12|), with N 
given by the general expression (g), the Euler-Lagrange 
equations give the amplitude a and width (3 



a 2 = (A + j D//3 2 )(2 + 2ct 2 //? 2 ) d / 2 , 



(3 2 = [4 - D + V(4-£) 2 + 16Aa 2 ]/(2A), 



(13) 



In Fig. H we show the power P s ~(tt /2) D / 2 a 2 j3 D and 
Hamiltonian of the stationary solutions in 2D. The 
dashed lines give the results of the weakly nonlocal ap- 
proximation with N given by Eq. (0), from which a 2 = 
4A and 1 = 2/A + 2cr 2 is found, resulting in the power 



P s = 4tt(1 + a 2 A), 



(14) 



where 4ir is the (A-independent) power of the Gaussian 
approximation to the soliton solution of the standard 2D 
NLS equation, recovered in the local limit ct=0. 
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FIG. 2. 2D variational results with Gaussian response and 
trial function. Left: Soliton power (solid) versus eigenvalue 
A for different degrees of nonlocality, <r=0, 0.2, 0.4, and 0.6. 
Dashed lines show the weakly nonlocal approximation. Right: 
Corresponding Hamiltonian versus power diagrams. 

In the 2D NLS equation the collapse is critical and 
the stationary solutions are "marginally stable" with 
dP s /d\=0 |l],py§. Typically perturbations act against 
the self-focusing, with several effects, such as non- 
paraxiality and saturability, completely eliminating col- 
lapse [||. This is also the case with nonlocality, as ev- 
idenced from Fig. [2] and the simplified expression (pd|), 
which shows that any finite width of the response func- 
tion (non-zero value of <j) implies that dP s /d\ becomes 
positive definite. According to the (necessary) Vakhitov- 
Kolokolov (VK) criterion |26j the soliton solutions there- 
fore (possibly) become linearly stable. 

For small A the soliton width (3 decreases as 1/A. Thus 
the assumption of weak nonlocality, i.e. that the soliton 
is much wider than the response function, applies only to 
sufficiently small values of A satisfying Act 2 <C 1, which 
is also clearly seen from Fig. |[ The accuracy of the 
assumption of weak nonlocality is further discussed in 
Ref. [^9| in terms of modulational instability. 

The 3D case shown in Fig. || is more interesting, be- 
cause the nonlinearity is much stronger than in 2D. The 
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collapse in the local 3D NLS equation is so-called super- 
critical Jl|,§,||. Again the soliton width (3 decreases as 
1/A, so a threshold width should exist, below which the 
nonlocality is not strong enough to stabilize the soliton. 
This is exactly what is observed in Fig. |^: For A < A th the 
solitons are still linearly unstable with dP s /d\ < 0, but 
above threshold the nonlocality is strong enough to bend 
the curve and make dP s /dX > 0, i.e., the solitons be- 
come linearly stable according to the VK criterion. From 
the definition dP s (X th ) j d\=0 the variational results give 
A th = l/(2(7 2 ), corresponding to a threshold in the soli- 
ton power (dashed curve in Fig. |^) and width 

P s th = (5^) 3/2 5ct/4, /3 th = 2a, (15) 

which are both proportional to the degree of nonlocal- 
ity a. Thus, sufficiently broad and high-power solitons 
are stable. In the Hamiltonian versus power diagram 
in Fig. H the lower (upper) branches correspond to stable 
(unstable) solutions while the threshold is represented by 
the cusp p0|| . This stable solution branch was recently 
found numerically in the context of BEC's with a nonlo- 
cal negative scattering potential p2]| . It corresponds to 
high density, self bound states of the condensate. 
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FIG. 3. 3D variational results with Gaussian response and 
trial function. Left: Soliton power (solid) versus eigenvalue 
A for different degrees of nonlocality, a—Q, 0.2, 0.4, and 0.6. 
Dashed lines show the threshold power ([l5|). Right: Corre- 
sponding Hamiltonian versus power diagrams. 

In conclusion we studied the properties of localized 
wave packets in nonlocal NLS equations. We proved 
rigorously that nonlocality of an arbitrary shape elimi- 
nates collapse in all physical dimensions. We also demon- 
strated that multidimensional soliton solutions of the 
NLS equation may be stabilized by the nonlocality. This 
opens a new and interesting discussion as to what is actu- 
ally observed in collapse experiments in nonlocal systems. 
It seems clear that it all comes down to oscillations be- 
tween opposite extreme states and how strong and rapid 
they are. Such oscillations were recently found to occur 
in BEC's through numerical and variational studies 

We acknowledge useful discussions with Yu.S. Kivshar, 
S.K. Turitsyn, and Yu.B. Gaididei. This work was sup- 



ported by the Danish Technical Research Council (Talent 
Grant No. 26-00-0355) and the Danish Natural Sciences 
Foundation (Grant No. 9903273). 



[1] J.J. Rasmussen, K. Rypdal, Phys. Script. 33, 481 (1986). 

[2] L. Berge, Phys. Rep. 303, 259 (1998). 

[3] Y.S. Kivshar, D.E. Pelinovsky, Phys. Rep. 331, 117 

(2000) . 

[4] A.Y. Wong, P.Y. Cheung, Phys. Rev. Lett. 52, 1222 
(1984). 

[5] E. Garmire, R.Y. Chiao, C.H. Townes, Phys. Rev. Lett. 

16, 347 (1966). 
[6] C.A. Sackett et al, Phys. Rev. Lett. 82, 876 (1999). 
[7] P.P. Banerjee, A. Korpel, K.E. Lonngren, Phys. Fluids 

26, 2393 (1983). 
[8] S.L. Shapiro, S.A. Teukolsky, Black Holes, White Dwars, 
and Neutron Stars - The Physics of Compact Objects 
(John Wiley & Sons, New York, 1983). M. Houbiers, H.T. 
Stoof, Phys. Rev. A 54, 5055 (1996). 
[9] C. Sulem, P. Sulem, The Nonlinear Schrddinger Equa- 
tion, Springer- Verlag, Berlin (1999). 
[10] V.E. Zakharov, in Handbook of Plasma Physics Vol. 2. 
Basic Plasma Physics II, eds. Galeev and Sudan (North- 
Holland, 1984). 
[11] F. Dalfovo et al., Rev. Mod. Phys. 71, 463 (1999). 
[12] O. Bang et al., Phys. Rev. E 49, 4627 (1994); Applicaple 

Analysis 57, 3 (1994). 
[13] O. Bang, D. Edmundson, W. Krolikowski, Phys. Rev. 

Lett. 83, 5479 (1999). 
[14] A.G. Litvak, JETP Lett. 4, 230 (1966). 
[15] D. Suter, T. Blasberg, Phys. Rev. A 48, 4583 (1993). 
[16] A.G. Litvak et al., Sov. J. Plasma Phys., 1, 31 (1975). 
[17] T.A. Davydova, A.I. Fishchuk, Ukr. J. Phys. 40, 487 
(1995). 

[18] H.L. Pecseli, J.J. Rasmussen, Plasma Phys. 22, 421 
(1980). 

[19] L. Berge, A. Couairon, Phys. Plasmas 7, 210 (2000). 
[20] K. Goral, K. Rzazewski, T. Pfau, Phys. Rev. A 61 

051601R (2000). 
[21] V.M. Perez-Garcia, V.V. Konotop, J.J. Garcia-Ripoll, 

Phys. Rev. E 62, 4300 (2000). 
[22] A. Parola, L. Salasnich, L. Reatto, Phys. Rev. A 57, 

R3180 (1998). 
[23] S.K. Turitsyn, Teor. Mat. Fiz. 64, 226 (1985). 
[24] D.V. Fedorov, A.S. Jensen, Phys. Rev. A 63, 063608 

(2001) . 

[25] W. Krolikowski, O. Bang, Phys. Rev. E 63, 016610 
(2001). 

[26] X. Wang et al, Phys. Rev. A 37, 3557 (1988). 

[27] A.W. Snyder, D.J. Mitchell, Science 276, 1538 (1997). 

[28] N.G. Vakhitov, A. A. Kolokolov, Radiophys. and Quan- 
tum Electronics 16, 783 (1975); 

[29] W. Krolikowski et al, Phys. Rev. E 64, 016612 (2001). 

[30] N. Akhmediev, A. Ankiewicz, R. Grimshaw, Phys. Rev. 
E 59, 6088 (1999). 



4 



